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Abstract. We provide a brief review of some of the major research results 
arising from the method of the Inverse Scattering Transform. 



1. Introduction 



I will give a brief review of several items in the Legacy of the Inverse Scattering 
Transform. In no way is this to be a complete review, since the Legacy has become 
so vast. However, I will treat those items with which I am most familiar, and try 
to detail their significance and importance. 

There is no doubt that the most important contribution was the famous clas- 
sical Gardner, Greene, Kruskal and Miura (GGKM) work Q of 1967 on the KdV 
equation. This was the starting point. They had found a very strange and new 
method for solving the initial value problem of a nonlinear evolution equation, the 
KdV. At that time, and even for several years later, this strange new method was 
considered to be only a novelty, since it would only work for that one equation, the 
KdV. Shortly thereafter, as a prelude to what was to follow, Peter Lax (2j showed 
that if given an appropriate linear operator, L, dependent on a potential, u(x), then 
one could always construct an infinite sequence of evolution operators, B, each of 
which would satisfy 



This sequence of evolution operators could be generated by simply increasing the 
order of the spatial differentials contained in B. Then from (1.1) one would obtain 
additional nonlinear evolution equations, each of the form 



where K was some (nonlinear) operator. All these additional higher order evolution 
equations would be solvable by this same technique. This collection is now known 
as the KdV hierarchy. 
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BL - LB — d t L. 



(1.2) 



d t u = K(u) 
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If we consider the eigenvalue problem for L, 
(1.3a) Li) = -Xtp, 

where A is the eigenvalue, and append to it the condition 
(1.3b) dti> = Bip, 

then it is easy to see that (1.1) is simply the integrability condition for (1.3). 
Furthermore, as Lax pointed out for the KdV case, (1.1)-(1.3) also implies that the 
eigenvalues, A, in (1.3a) would be stationary, 

(1.4) d t x = 0, 

a relation that would occur time-and-time again as the study of integrable equations 
would expand in the decades to follow. 

It was also about this time that the term "radiation" was introduced. We 
haven't said anything yet about solitons or solitary waves, but more will be said 
later. For now, let us note that a remarkable feature of the GGKM method of so- 
lution was the appearance of fully nonlinear solitary wave solutions, called solitons. 
The other part of the solution has been called "radiation" , and is essentially linear- 
like in its behavior. The asymptotics (long-time behavior) of the total solution are 
generally that the radiation does disperse away, leaving the solitons traveling in a 
sea of decaying radiation. 

As to nomenclature, we shall refer to (L + X)tp as the eigenvalue problem, ip as 
the eigenf unctions, and B as the evolution operation. The pair [L + A, B] is known 
as the "Lax pair" . For the KdV equation, the Lax pair is 

(1.5a) (L + \)i/j= (d 2 x + u + \)i) = 

(1.5b) d t ip = hl>={a- Adl - 6ud x - 3u x )i> 

where a is an arbitrary constant and the integrability condition is the KdV equation: 
(1.5c) d t u + d x u + 6ud x u = 0. 

We note here, given L and B, it follows that one can then obtain K(u). How- 
ever, an important problem is given K(u), construct L and B. The solution of 
this inverse problem is still an area of active research. One method that sometimes 
works for this is called "Painlcve Analysis" . For a description of this aspect of the 
Legacy, the reader is referred to Choudhury's article in this same issue. 

It was not until 1971, that the next physically significant integrable system 
was uncovered by V. Zakharov and A.B. Shabat (ZS) Q, which was the focusing 
Nonlinear Schrodinger Equation (NLS) 

(1-6) iq t = ~q X x -2q*q 2 . 

This equation required a different eigenvalue problem, 
(1.7a) vi x + i(vi = qv 2 

(1.7b) v 2x - Kv2 = rvi 

where (v\, v 2 ) is the eigenvector, £ is the eigenfunction, and q and r are the "poten- 
tials". For the focusing NLS case, one has r = — q* and r = +q* for the defocusing 
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case. The time evolution operator, B, is given by 

(1.8a) id t Vx = -i(4( 2 + 2qr)v 1 + (A(q + 2iq x )v 2 

(1.8b) id t v 2 = (4Cr - 2ir w )v l + i(4( 2 + 2qr)v 2 . 

Exactly as was shown by Lax [§] for the KdV, one also has a hierarchy here, which 
can be obtained by generalizing (1.8) to higher orders in £. In 1972, Wadati 
found the next member of this hierarchy, the "modified KdV" (mKdV) 

(1.9) q t + q X xx + Gq 2 q x = 

which was also integrable. Its eigenvalue problem was again the ZS eigenvalue 
problem, (1.7), but where now r = —q, and q real. Also, (1.8) had to be generalized 
to be cubic in £. 

By this time, it was becoming apparent to many researchers, that this strange 
method found by GGKM was not simply a novelty. Rather, there was some- 
thing very significant underlying all of this. This became even more obvious when 
Ablowitz, Kaup, Newell and Segur (AKNS) presented a method of solution of both 
the Goursat and Cauchy initial value problems of the sine-Gordon equation |§. 
This was also based on the ZS eigenvalue problem, but with a very different form 
for the B in (1.8): it was now inversely proportional to the spectral parameter, 
(. The sine-Gordon equation was well known at that time. It had a long history, 
first occuring in 1853 in differential geometry, and was the first equation for which 
Backlund tranformations and N-solitons solutions were found. It was known in 
solid state physics in the 1930's, and in 1965 had found applications in optics. 

The 1ST solution of the sine-Gordon equation was shortly followed by another 
letter Q pointing out how one could generate a large number of integrable equa- 
tions, each of which were physically significant and important. With one general 
approach, AKNS were able to reproduce all the Lax pairs found up to that time, 
and were able to connect the form of the dependence of B on £ to the linear dis- 
persion relation, uj(k). (The linear dispersion relation relates how the frequency, 
u), depends on the wave vector, k, in the linear limit, where plane waves, e l '~ kx ~ u ' t ', 
are the natural solutions.) In 1974, they published their classic AKNS paper Q, 
wherein they described in detail this new method of solution, calling it the method 
of the Inverse Scattering Transform (1ST). One of the major points of this classic 
was that the 1ST could be viewed as a nonlinear extension of the method of the 
Fourier Transform. 

This was also the start of the explosion in research on solitons and integrable 
systems, because unbeknownst to most westerners, Faddeev, Zakharov and their 
students were all very busy in the same direction. In the next few years, many 
important papers were to be published on the 1ST and related issues. 

2. The Legacy 

Beginning in 1974, it becomes difficult to try to detail all the results. Never- 
theless, we will now discuss in general terms, the legacy which followed from this. 
In the following, we will list the general areas of the legacy, and briefly describe the 
importance and the major contributions made to each one. 

2.1. Method of Solution — the 1ST. Above all, the 1ST is a method of so- 
lution for integrable nonlinear equations. It was the pioneering work of GGKM Q, 
ZS [U and AKNS Q] which made the most significant impact and set the tone which 
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followed. Consider what the 1ST allows one to do. One may take any reasonable 
initial data, and by means of the direct scattering transform (the eigenvalue prob- 
lem) , transform the initial data into "scattering data" . For the KdV, the scattering 
data consists of a reflection coefficient [p(k); — oo < k < oo], bound state eigenvalues 
[kj; j = 1,2, ... ,N] and bound state normalization coefficients [Cf,j = 1, 2, . . . , N] 
where N is the total number of bound states, usually finite. Now, Lax Q showed 
that the eigenvalues, Kj, would be independent of time, due to (f .1). GGKM showed 
that if u(x,t) evolved according to the KdV equation, the reflection coefficient, p, 
and the normalization coefficients, Cj, would evolve according to 

(2.1a) d tP (k;t) = p(k;0)e Stk3t , 

(2.1b) d t Cj(t) = C J (0)e- 8K K 

Thus it becomes a very simple matter to determine the scattering data at any later 
time. 

Next, one used the solution of the inverse scattering problem, which is the core 
of this method of solution, to reconstruct the potential(s). One transforms (with 
the 1ST) from the scattering data, at time t, back to the potential(s), at time t. 
For the KdV, the necessary steps are to first construct 

1 roc N 

(2.2) F{z;t) = — P {k;t)e ila dk + Y t Ci{t)e-^', 
then one solves the linear integral equation 

(2.3) K(x,y) + F{x + y)+ / K(x, s)F(s + y)ds = 0, 

J X 

for K(x,y;t). Lastly u(x,t) is constructed from 

,„ . . „dK(x,x:t) 

(2.4) u(x,t) = -2 V-^. 

ax 

All integrable systems are solved by an 1ST of the above format, although there can 
be a wide variation in the form of the formats. Some eigenvalue problems are higher 
order and/or even multidimensional. But there is always some scattering problem 
which maps the potential(s) into a set of scattering data (the Direct Scattering 
Transform). There is always some evolution of the scattering data, as in (2.1). 
There is always an Inverse Scattering Transform that allows one to map from the 
scattering data back to the potential(s), as in (2.2)-(2.4). 

One may describe this method of solution as a method for solving the initial 
value problem of a nonlinear equation by using only linear techniques. Further- 
more, one may also say that one solves these nonlinear problems by not solving the 
nonlinear problem. Instead one solves two related linear problems. 

To understand these comments, consider 

(2.5) L(u,C)-V = 0, 

where L is the eigenvalue problem in the Lax pair, as in (1.1)-(1.3) and £ is a spectral 
parameter. Let the second component of the Lax pair be B{u, u x , . . . , d x , . . . ), 
where 



(2.6) 



B - V = d t V. 
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Now note that (2.5) determines the ^-dependence of V(x,t;Q, whereas (2.6) de- 
termines the i-dependence. Thus one function, V(x,t;Q), is being determined by 
two equations. In general, such would overdetermine V, and in order for a mutual 
solution to exist, certain consistency conditions must be satisfied. This condition 
is the single integrability condition, (1.2). When (1.2) is satisfied, then a mutual 
solution exists for (2.3)-(2.6). 

But if this integrability condition is nothing more than the nonlinear evolution 
equation, it therefore follows that if, by some means, we can construct a single- 
valued solution for V(x,t;Q, for some u(x,t), which satisfies each component of 
the Lax pair, then it follows that u(x,t) must satisfy the nonlinear evolution equa- 
tion. So, one could say that the entirety of the method of the 1ST is based on not 
solving integrable nonlinear equations (at least, not directly). Instead, we solve 
them indirectly, exactly by the same format used in any transform method. This 
is as follows. We satisfy (2.5) by mapping u(x, i) into scattering data and con- 
structing the eigenfunctions. We can do this for certain classes of potentials (i.e. 
L\f)L,2). Then we always will have reflection coefficient (s), as a function of the 
spectral parameter, £, and certain bound state data (eigenvalues and normalization 
coefficients). This is a linear problem. We satisfy (2.6) by requiring the scatter- 
ing data to evolve appropriately (as in (2.1) for the KdV). This is usually easy to 
do, since scattering data is typically defined for x — ► ±oo, where potentials nor- 
mally approach specified values (usually zero) in that limit. So one really only needs 
B(x — ► ±oo). This is also a linear problem. Therefore, by fixing the scattering data 
to evolve appropriately, we have effectively forced u(x, t) to evolve by the nonlinear 
evolution equation. To solve for u(x,t), we need to solve the inverse scattering 
problem, which is also a linear problem. The construction of the kernel(s), F(z), 
as in (2.2), is another linear problem. Equation (2.3) is a linear integral equation 
for K(x,y). Then u is reconstructed as in (2.4) by a linear operation. Thus with 
only linear techniques, we are able to solve these nonlinear evolution equations. 

2.2. Soliton Solutions. One of the unique features of the 1ST is that it 
allows one to construct an infinity of exact nonlinear solutions, called the TV-soliton 
solutions. These are also called reflcctionless potentials because the scattering data 
consists of only bound state scattering data, with all reflection coefficients set equal 
to zero. In the case of the KdV equation, when the reflection coefficients vanish, 
the function F{x + y;t) in Eq. (2.3), then separates into a finite sum of products 
of known functions of x and y, allowing one to obtain K(x, y), also as a finite sum 
of known functions. The same is true for all other cases integrable by the 1ST. 

The value of these solutions is tremendous, since they allow one to study and 
obtain exact results for these systems. The most important solution in this class is 
always the one-soliton solution, since it is the basic building block of these solutions, 
and also of any interactions between solitons and between solitons and radiation. 
The 1-soliton solution of the KdV hierarchy is 



(2.7) 



2 V 2 



u = 



cosh 2 {77 [:r — xo(t)]} 



and that of the NLS hierarchy is 



(2.8) 




Q = 



cosh{2?7[x — £()(*)]} 
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where the form of xo(t) and <fio(t) depends on the member of the hierarchy. The 2- 
and 3-solitons solutions demonstrate soliton interactions and collisions. Typically 
what happens in any soliton collision is that after the collision, the solitons separate 
out according to their individual velocities, completely unscathed, except for a 
possible phase shift in their positions and/or phase. There are two well known 
exceptions to this. In the soliton decay case of the 3WRI [Q, an initial soliton 
in the high frequency wave can decay into its two daughter waves, transferring 
its identity to them. In the vector NLS [0], similar inelastic type collisions occur 
whereby a soliton in one mode (polarization), can flip to another combination of 
modes. Bound states can occur in the sine-Gordon field, where one can have stable 
2-soliton bound states, called "breathers" [^|, which are localized oscillations of the 
sine-Gordon field. 



2.3. Hierarchies. Given any eigenvalue problem, Lax j2[ had noted that one 
could always take the eigenvalue problem (the first component of the Lax Pair), and 
by simply extending the order of the evolution operator, B, one could generate an- 
other nonlinear integrable equation. All these nonlinear integrable equations which 
have a common eigenvalue problem in the Lax Pair, but different evolution oper- 
ators, B, is called a "hierarchy". Thus, there is a hierarchy for every one of these 
eigenvalue problems. For the Schrodinger equation, the most important members 
are the KdV equation, a 5th order KdV equation 
ton" equation 



11 



10 1 and one-dimensional "cavi- 



(the analogy of the sine-Gordon equation for the Schrodinger 
eigenvalue problem). In the ZS problem, we have a "workhorse" as far as physical 
equations are concerned. If we allow r in (1.7) to be in general independent from q, 
then in addition to the NLS JH, the hierarchy contains the modified KdV Q], the 
sine-Gordon JU , the sinh-Gordon equation H , coherent pulse propagation and self- 



induced transparency (SIT) |12|, stimulated Raman scattering (SRS) 13 ] , and the 
defocusing NLS [14]. The hierarchy containing the three- wave resonant interaction 
(3WRI) includes all three forms of this interaction (explosive, soliton decay, and 
stimulated backscatter (SBS)) 



15 



16], as well as the Manakov vector NLS 



2.4. Inverse Scattering. Another aspect of the legacy is the wide variety of 
inverse scattering problems solved. In 1967, the 1ST of the Schrodinger equation 
had just recently been obtained [17, 18], and it was only in 1972 that the 1ST of 
the ZS eigenvalue problem had been solved Since that time, there has been a 
multitude of other and even more complex scattering problems solved. 

The next one was the solution of the t hird -order eigenvalue problem for the 
three wave resonant interaction (3WRI) [15, 16 , which interaction we shall return 



to later. This one differed from the ZS significantly only in the additional order 
of the problem. An important subcase of this was the inverse scattering solution 
for the eigenvalue problem for the vector nonlinear Schrodinger (VNLS) equation, 
solved by Manakov Q] . This VNLS equation is a very important and key equation 
for several studies of nonlinear optical pulses propagating in optical fibers |19, 20 1. 



The generalization of the order 3 problem to order n has been done by Gerdjikov 
and Kulish 



21 



More complex forms have also been done. The first one of these was the eigen- 
value problem for the sine-Gordon equation in laboratory coordinates 



22, 



23 



Here one has the spectral parameter distributed among various potential terms, 
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and is 



(2.9a) 
(2.9b) 



1 i 

M I TTs - — " ) Vi 

1 i 

l '- J 1 2^ ~ 8C C ° SU ' Vl 



i . 1 , > 

— sin u - ~(u x + it*) 
oC, 4 

* . 1 / \ 

— sin u - -iu x + u t ) 



V-2 



Vl. 



One notes that due to this form, one cannot easily express this equation in the 
standard form as L ■ V = XV where L is a nondegenerate differential operator, V 
is the eigenvector and A is the eigenvalue. However, in spite of this, one still can 
solve the direct and inverse scattering problems for such systems. Further examples 
of such eig envalue problems are the eig env alue problems for the massive Thirring 
model [241, and the derivative NLSL |25|. A more standard form is the cubic 



generalization of the Schrodinger equation [1C, 26 1, 



(2.10) 



ipxxx + 6Q4>x + 6i?^ = AV> 



which appeared as the eigenvalue problem for the Boussinesq, Sawada-Kotera equa- 
tion and the Kaup-Kuperschmidt equation. There is also the inverse scattering for 
multidimensional problems, such as the 3D form of the 3WRI [27|, the KPI and II 
equations |28], as well as the DSI and II equations [29|. Several aspects of these 



are still of current research interest. 

2.5. Perturbations and Closure. Once one has an exact method for the 
solution of a system, it then becomes possible to develop perturbation methods, to 
study nearby systems. This work was first done in 1976 for the ZS eigenvalue prob- 
lem for the one-soliton solution 



30 



with a general summary of the perturbation 

method being given in 1978 

Key to this, is the concept called "closure" or "completeness" , which itself arises 
from the one-to-one nature of the direct scattering transform and the 1ST 



32 



What this simply means is that given any potential in the appropriate class, there 
exist a unique set of scattering data that can be associated with it by the eigen- 
value problem of the Lax pair, and vice versa for the 1ST. Thus for any potential 
whose evolution is slightly perturbed away from its integrable value, by the direct 
scattering transform, it will be mapped into some other scattering data near the 
initial integrable scattering data. The perturbation problem is then to determine 
how this scattering data in the perturbed case evolves in time. Once its evolution 
is known, then by the 1ST, one may map back to the potentials and then obtain 
their evolution. 

One solves this by relating variations in the potentials to variations in the 
scattering data, with the transformation from the former to the latter being ac- 
complished by the so-called "squared eigenfunctions" . Then with these squared 
eigenfunctions, one may obtain the evolution of the scattering data under a given 
perturbation. 

Under perturbations, one no longer has the simple evolution of the scattering 
data as in (2.1). Rather, one has a slow mixing of the various elements of the 
scattering data: solitons will decay and/or be pumped, transforming some of their 
energy into radiation and/or absorbing energy from the perturbations. In addition, 
radiation modes will similarly grow and/or decay. Many aspects of this have been 
covered in the review of soliton perturbations done by Kivshar and Malomed ]33| 
as well as in Ref . [31]. 
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However, whenever the eigenvalue problem has a singular structure, one has 
difficulties. The classical example is the perturbation theory of the KdV, where the 
Schrodinger eigenvalue problem is singular at k — 0. What happens is that for a 
pure one-soliton solution, p(k) = at k — 0. But if there is just even the smallest 
amount of radiation present, we have p(k = 0) = — 1. Now, the width of this 
region in fc-space can be very small (proportional to the radiation density) , so this 
region could be vanishingly small and one might have to look very close to k = 
in order to see it. However, due to this singular behavior, any perturbations of the 



KdV soliton which will create any radiation, will always have secular terms 34 . 



Physically what is occurring is that a shelf or a depression is forming due to the 



perturbation [31], and this action does generate a finite shift in the scattering data 



as the perturbation vanishes. 

2.6. General Integrable Evolution Equations. Given a linear dispersion 
relation, what are the possible integrable evolution equations for that system? Well, 
for any integrable system solvable by the ZS 1ST or the Schrodinger 1ST, that answer 



has been given by Newell and Kaup [35]. They showed that the most general system 
is any of the AKNS polynomial forms Q , coupled with a generalized SIT system. 
Using the properties of the ZS squared eigenstates and their closure, they were 
able to construct the most general evolution equation, given the linear dispersion 
relation. In other words, the dispersion relation of the linear theory determines the 
nonlinear terms, given the eigenvalue problem. The same could be done for any 
other hierarchy. 

Now, a burning question has always been "Is it ever possible for an exactly inte- 
grable system to have evolving eigenvalues (i.e. for soliton amplitudes to evolve)?" 
In general, the answer to this is "No" . However, Kaup and Newell did find ex- 
ceptions. Such equations can indeed be constructed, but potential applications for 
them seem to be remote. 

2.7. Optical Systems. One question which arose naturally after a few years 
is "Why are so many of the physical integrable systems related to nonlinear optical 
systems?" As one goes down the list, one has the sine-Gordon equation (which 
applies to two-level atoms and is the sharp line limit of SIT), nonlinear Schrodinger 
(focusing and defocusing — the NLS is almost always the weakly nonlinear limit of 



any almost monochromatic envelope [36 1), self- induced transparency (SIT — also 



more generall y r eferred to as two- level coherent propagation) , thr ee-level coherent 



propagation 37 1, three- wave resonant interactions (3WPJ) |15_L [L6L m, second 



41 



harmonic generation (SHG) |38|1 , the three-dimensional form of the 3WRI [|39| , 40, 



stimulated Raman scattering (SRS) |p.3| , |42j , |43|], two photon propagation 



(TPP) [£4 £5j, and degenerate TPP (DTPP) |47j . This is not meant to be 
a comprehensive list, but it does include the major integrable nonlinear optical 
systems. 

Of these, the 3WRI was the first system to demonstrate a major departure 
from the accepted and expected soliton behavior. The first deviation was that the 
radiation would never asymptotically vanish. Also after any and all collisions, nor 
would the radiation separate out from the solitons. The reason for this is that 
the 3WRI has no dispersion, and therefore solitons will never separate from the 
radiation present. They each have the same velocity. One could now ask what is 
the significance of the 3WRI solitons, when they differ so much from KdV or NLS 
solitons? To answer that the best we can say is that in the 3WRI, any solitons 
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seem to simply represent a "packet" of something. It is a unit which cannot be 
broken up. Although there are certain exchange rules for the exchanging of solitons 



between the three envelopes [16, |8|], nevertheless solitons in the 3WRI do seem to 
be some robust and coherent part of the envelope. On the other hand, the radiation 
component of any envelope has no such finite "packet" size, but rather can always 
be subdivided and redistributed among the three envelopes, subject only to the 
Conservation of Action laws Q . 

It was the SRS system that first brought to the forefront certain ambiguities 
with the 1ST on finite intervals. The first solution of the finite interval case, by 



using an infinite interval 1ST, was given in [45], where the general 1ST for SRS 



was developed, and features of the solution were discussed. Numerics of SRS have 



been studied by Hilfer and Menyuk 48 1 , the asymptotical form of the solution was 
described by Kaup Et|, and was later solved as a Riemann-Hilbert problem by 
Fokas and Menyuk | 43[ . In the meantime, by analytically solving a model initial 
value problem, Menyuk and Kaup essentially found that for all these integrable 
nonlinear optical problems (except NLS and 3WRI), one could just as easily describe 
the solution as "an infinity of solitons with no radiation", as well as by the usual 
description of "a finite number of solitons in a sea of radiation" . Briefly, why this 
could happen is basically the same reason as why a Fourier transform on a finite 
interval has a variety of forms. There one could use either a cosine series, or a sine 
series, or the exponential series to represent the function. For a given function, the 
coefficients in each of these series is quite different. Another way to look at this, is 
that on a finite interval, there is an infinite number of ways to take and combine 
plane waves to reconstruct some function inside the finite interval . Similarly for 
the 1ST, on a finite interval, there is no unique form for the scattering data. 

Pursuing this further, consider the case where solitons are forbidden, as in the 
ZS r = +q* case. Now, what is going to happen in this case where no solitons 
are allowed on the infinite interval, if we try to represent the solution on the finite 
interval with no reflection coefficients, and only with solitons? What happens to 
this system on a finite interval? Well, we again find something surprising. Taking 
SHG as an example case of these defocusing systems, we find that the solution 
of a simple initial value problem can also be given in terms of no radiation, but 
now (since regular solitons are forbidden) an infinity of what is called "virtual 
solitons" |50fl . These are indeed interesting objects. 

The possibility of their existence was noted way back in 1974 Q. However, no 
known application was then known for them. In the ZS 1ST, the typical soliton for 
the focusing case (r — —q*) is of the form 

2fj e - 2i i( x ~ x o) e -*0o 

^ ^ cosh[2?7(a; — .To)] 

Now, if one simply assumes that the defocusing case (r — +q*) has one bound state 
in the scattering data, one obtains the "virtual soliton" solution 

2'/ 



(2-12) q = 



ih[2r)(x - x )] 



which clearly will always be singular on any infinite interval. However, on a finite or 
a semi-infinite interval, as long as Xq lies outside the physical region, the solution 
is perfectly valid. Thus a virtual soliton is the singular, ZS, r = +<?*, soliton. 
(For certain technical reasons, they also were first found in the lower half complex 
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plane [50], which is another reason for the prefix "virtual". For more details on 
this, see Steudel's contribution in this same series.) 

As to my first question as to why so many nonlinear optical systems are in- 
tegrable, perhaps the principle reason for this is the vast orders of magnitude dif- 
ference between the speed of light and other velocities in these physical systems 
(which are typically acoustic velocities, and/or electronic drift velocities). Due to 
this disparity in velocities, one can then expect a multiple-scale expansion to give 
quite good results. Also one would then expect to see the higher order terms scaling 
like some power of this velocity ratio, and therefore rapidly vanishing. 

2.8. Benjamin-Ono Equation. Although we could perhaps have included 
the Benjamin-Ono (BO) equation in the section on eigenvalue problems, it is unique 
enough to justify some additional comments. First it is a one dimension eigenvalue 



problem but it is also a nonlocal eigenvalue problem |51|. Actually it can be 
formulated as an electrostatics problem (Poisson's equation) in two dimensions, 
since it can be stated in terms of Hilbert transforms and their analytical properties. 



This is probably the simplest explanation as to why Fokas and Ablowitz [51] could 
term the BO equation as a "pivotal equation for multidimensional problems" . It 
indeed does contain this multidimensional flavor. It also should be noted that there 



are now two versions of its 1 ST [51, 52|. This is not entirely surprising, since the 
multidimensional 3D-3WRI 41 1 also contains a multitude of different forms of the 
1ST. The multidimensional flavor of the BO equation stands out even more when 
one considers its closure property, its perturbation theory, and its Hamiltonian 



structure 53 



2.9. Reduction Problems. It seems that almost any integrable system can 



be found in the Yang-Mills field, if one knows how to find the right reduction 54 1. 
However, there is another quite useful direction that one can take in reducing 
integrable systems. Let us first note that if one simplifies, or reduces the number 
of degrees of freedom of an integrable system by some set of constraints, consistent 
with the integrable flows, then the reduced system will also be integrable. As one 
example, consider the propagation of N solitons in an optical fiber. It is very 
important to maintain the spacings of these N solitons over long distances. So 
a key question is the stability of such an arrangement, which is simply an N- 
soliton solution. However, to try to analyze a 100-soliton state, where every soliton 
has approximately the same amplitude, is not an easy thing to analyze, since the 
transmission coefficient, a, has a zero of order 100 at the pole. Another approach 
for N larger than 2 or 3 is clearly needed. One way to do this is to approximate 
the system as a lattice, since in general, one expects the solitons to have an almost 
equally spacing (or being absent if representing a zero) and of equal amplitude. 
Then it turns out that when the A-soliton system is reduced to a lattice system, 
it reduces to a Complex Toda Lattice of N points, which is integrable. One now 
can study the stability problem of the reduced system, the Complex Toda Lattice, 



and transfer the results to the A-soliton optical pulse [55]. From this, one obtains 
general criteria about how the phase of each successive soliton should be adjusted, 
to maximize the stability of the pulse train, and etc. 

2.10. The Fokas Method. Another new aspect of the 1ST has been recently 



developed by Prof. Fokas [56]. There are complex practical problems requiring 



the solution of 2-dimensional, linear, constant coefficient, PDEs but with complex 
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boundaries, such as wedges and polygons. Before his work, there was no general 
analytical technique for constructing solutions with such complicated boundaries. 
However, it is a simple matter to create a Lax pair for almost any linear PDE with 



constant coefficients [57|. Once this is done, one then can approach the solution 
of these linear two-dimensional, boundary value problems from the 1ST point of 
view, whereby one solves (or satisfies) both components of the Lax pair, thereby 
also satisfying the integrability condition, which is just the PDE to be solved. The 
method can also be extended to integrable nonlinear PDEs and evolution equations, 
however the solution then frequently requires the solution of a Riemann-Hilbert 
problem for the reflection coefficients. 

2.11. Unsolved Problems. Now we will discuss what the hard problems are. 
They are hard because they haven't been solved. This is not a complete list, but 
is the start of a list of problems in need of solutions. 

Although we hae a solution for the 1ST of the cubic eigenvalue problem, (2.20), 
there still is no equivalent of the GLM equations for this system. What we do have 
is the solution of the Riemann-Hilbert problem for the eigenf unctions. What is next 
needed is a representation of the eigenfunctions in terms of some transformation 
kernels (like the K(x, y; t) for the KdV). Once these are known, then the equivalent 
GLM equations will follow. Some progress in this direction has recently been made 
by A. Parker [§f]. 

There are still aspects of KPI and II that are of interest. Numerical simula- 



tions [ 5S ] provide valuable insight into the evolution of these equations. I would also 
say that there are probably still some unanswered questions about DSI and II. How- 
ever, the unfortunate thing about these two equations is that potential applications 
seem to be almost lacking, due to the scales involved. (Typical parameters for DS 
solitons in water require meter-like distances horizontally, but only centimeter-like 
water depths, and even smaller wave amplitudes.) 

Questions about perturbation theory and closure relations seem to have become 



quite well understood, based on the works of Gerdjikov, Ivanov and Kulish |60|, 



Beals and Coifman 61], and J. Yang [62|. They are also quite well understood for 



the 3D3WRI and DSI, since they both use the same eigenvalue problem. However, 
I know of no published work in this area. 

Of the optical problems there is still interest in SHG and DTPP, and the 
latter has even the direct scattering problem to be detailed, as well as the inverse 
scattering problem. 

I also strongly suspect that there are other integrable systems still to be found. 
Here one would have to apply Painleve analysis and see what will result. 
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